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We discuss the positive indications of a possible dark matter signal in direct detection experiments 
in terms of a mechanism of interaction between the dark matter particle and the nuclei occurring 
via the exchange of a light mediator, resulting in a long-range interaction. We analyze the annual 
modulation results observed by the DAMA and CoGeNT experiments and the observed excess of 
events of CRESST. In our analysis, we discuss the relevance of uncertainties related to the velocity 
distribution of galactic dark matter and to the channeling effect in Nal. We find that a long-range 
force is a viable mechanism, which can provide full agreement between the reconstructed dark matter 
properties from the various experimental data sets, especially for masses of the light mediator in 
the 10-30 MeV range and a light dark matter with a mass around 10 GeV. The relevant bounds on 
the light mediator mass and scattering cross section are then derived, should the annual modulation 
effects be due to this class of long-range forces. 



I. INTRODUCTION 

Dark-matter direct-detection experiments have been 
providing exciting results in terms of measured features 
which have the right properties to be potentiaUy ascribed 
to a dark matter (DM) signal. The typical effect of 
annual modulation of the recoil rate [T] has been put 
under deep scrutiny by the DAMA/Nal Collaboration 
starting more than a decade ago [2]. Annual modulation 
of the rate with viable DM interpretation was observed. 
The upgraded DAMA/LIBRA detector has confirmed [3], 
with much larger statistics, the annual modulation effect, 
reaching the unprecedented result of an evidence of 8.9 
a C.L. for the cumulative exposure [3]. The DAMA an- 
nual modulation effect has been shown to be compatible 
with a DM effect which, for the case of a coherent scat- 
tering, refer to a range of DM masses which spans from a 
few GeV up to a few hundred of GeVs and cross sections 
between 10""*^ cm^ to 10~^^ cm^ [2H3], and with some 
noticeable differences due to the galactic halo modeling 
016!. 

More recently, the CoGeNT experiment first reported 
an irreducible excess in their counting rate [7], which 
could be in principle ascribed to a DM signal. In the 
last months, the same experiment reported an additional 
analysis which shows that the time-series of their rate is 
actually compatible with an annual modulation effect [8] . 
The evidence of modulation in CoGeNT is at the level of 
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2.8 cr C.L.. 

The interesting feature is that the DAMA and Co- 
GeNT results appear to be compatible for relatively light 
DM particles, in the few GeV to tens of GeV mass range 
and coherent scattering cross section around 10""*^ cm^ 
to 10^**° cm^ [5] (as usual, the actual relevant range of 
masses and cross section depends on the assumptions on 
the galactic DM properties, namely the velocity distri- 
bution function and the local DM density [6]). Further 
relevant analyses can be found in Refs. [51-fT5] and Refs. 

[I4H26]. 

In this paper we discuss an alternative possibility, 
namely the case of a DM particle that scatters on the 
nucleus with long-range interactions, like those that are 
induced by a light mediator, such that the nature of the 
DM-nucleus cross interaction is not contact-like. An ex- 
ample of these kind of interactions is given by mirror 
photons in models of mirror dark matter [SJ . We 

analyze if and under what circumstances a long-range 
force can explain the positive hints of a signal in direct 
detection experiments and what kind of bounds could be 
derived by those evidences on the light mediator mass 
and scattering cross-section (or, alternatively, coupling 
constants). 

We discuss the impact of the light-mediator parame- 
ters (mostly its mass, which determines the level of devia- 
tion from the standard case of a contact-like scattering) 
on the reconstruction of the DM mass. We show that 
long-range forces mediated by a 10-30 MeV boson may 
provide compatibility between the different experimental 
direct-detection results. These results are discussed for 
some variation of the galactic halo models. 
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We concentrate our analysis mostly on the DAMA 
and CoGeNT results, since they are based on a spe- 
cific feature of the DM signal, namely the annual mod- 
ulation, which is hardly mimicked with the correct fea- 
tures (period, phase, energy range, size) by background 
sources. However we will discuss also the relevance of 
the CRESST irreducible excess of events [571 liH] , which 
is currently based on the total event rate (without resort- 
ing to time dependence) but nevertheless points toward 
an additional indication of a possible signal. 

CDMS and XENON experiments have recently re- 
ported a small number of events which pass all the se- 
lection cuts (2 events for CDMS [SI] and 6 events for 
XENON 100, reduced to 3 events after post-selection 
analysis UDDj still too few to be correlated to a signal. 
They therefore can provide upper bounds on the DM 
scattering cross-section. The actual response of these de- 
tector to the light DM case which is here under scrutiny 
has been critically analyzed and appears to be uncer- 
tain and model dependent for light DM [6TI - I63] . We will 
show the results obtained from CDMS and XENON in 
our analysis on long-range forces, but we will not enforce 
those bound. 



II. DIRECT DETECTION SIGNALS 

Direct detection relies on the direct scattering of DM 
particles off the nuclei of ordinary matter, the two main 
processes being elastic scattering: 

X+AA(A,Z)atrcst X+A/'(^,^)rocoih (1) 

and inelastic scattering: 



X + U{A, Z)at rest ^ X' + Af{A, Z)rccoil, 



(2) 



In Eqs. ([T]) and ^ x and x' are the dark matter parti- 
cles and its excited state and A, Z are the atomic mass 
and atomic number of nucleus A/", respectively. In the 
detector rest frame, a DM particle with velocity v and 
mass m^^, would produce a nuclear recoil of energy E^^. 
The minimal velocity providing a nuclear recoil energy 
Er is: 



/ niN Er 



XN 



1 



^J■xN 5 

mM Er 



(3) 



where 5 — m!^ — the mass splitting between x and x'- 
Elastic scattering occurs for 6 = 0, while S ^ implies 
inelastic scattering. 

In this paper we will consider only the case of elastic 
scattering, but we will extend the mechanism of interac- 
tion to the possibility that a light mediator may induce 
a long-range interaction, instead of the typical situation 



where the scattering cross-section is obtained through a 
contact interaction, like it is the case for the exchange 
of heavy bosons. Long range forces alter the detector 
response to DM interaction, as will be outlined in the 
remainder of the Section, and can make low-threshold 
detectors (like CoGeNT and DAMA) especially sensitive 
to the presence of these long-range forces. 



A. Generalization of the point-like cross section to 
long— range interactions 

Long range interactions can be described by means of 
a light (massless in the extreme limit) mediator (f) which, 
in the non-relativistic limit, suitable for the DM-nucleus 
scattering in DM direct detection, corresponds to the 
presence of a Yukawa potential, whose scale is determined 
by the mass of the mediator. A specific realization 
is offered by models of mirror dark matter and models 
where mirror photons possess a kinetic mixing with or- 
dinary photons [HI ETHiS] . In this case, mirror charged 
particle couples to ordinary nucleus with electric charge 
Ze (Z being the number of protons in the nucleus), with 
an effective coupling e Z'^dark (being Z'ljdark the coupling 
between DM and mirror photon, and with e parametriz- 
ing the kinetic mixing between mirror and ordinary pho- 
ton). The radial dependent Yukawa potential of the in- 
teraction can be cast in the form: 



y(r) = (asMttdark) 



xtZZ' 



(4) 



where asM = e^/(47r) is the electromagnetic fine struc- 
ture constant and c^dark — ffdark/C^"")- III a more gen- 
eral framework, one replaces [(asMadark)^e ^ with 
the relevant coupling factors between DM and the nu- 
cleus, which may be just on protons, or on protons and 
neutrons with suitable strengths determined by the spe- 
cific model. For definiteness, we consider here the case 
motivated by mirror photons, which implies interactions 
with protons only. 

From the potential in Eq. Q, one obtains the differ- 
ential cross section: 



da(y, Er) 



2r7ijvA 



> 2 „2 



F^{Er). 



(5) 



where = 2miqER is the square of the momentum 
transferred in the interaction, v is the speed of the DM 
particle, F{Er) denotes the nuclear form factor which 
takes into account the finite dimension of the nucleus 
and: 
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A 



27rQ!sMadarke^^^^'^ 



niN 



(6) 



The differential cross section of Eq. ^ exhibits two Hm- 
its. The point-like limit of the interaction occurs when 
the mass of the mediator is much larger than the trans- 
ferred momentum, i.e. when ^ q. By "point-like" in 
this context we mean that the mechanism of interactions 
is realized through a contact interaction. In this regime 
the differential cross section reduces to the standard case: 



da{v,ER) mjv 1 2 



dER 



2^2 „2 



Z a^^ F [Ei- 



(7) 



where /i^p is the DM-nucleon reduced mass and the total 
point-like cross section per nucleon for such coupling is 
given by: 



le-Tre^Z'^asMadark 



(8) 



Eqs. (TjSl can be generalized to the case of scattering 
off both protons and nucleons with a change such that 
— >• [Z + fn/ fp{A — Z)]'^, where A is the mass number 
of the nucleus and /p_„ are factors which differentiate the 
coupling on protons and neutrons. The case we are con- 
sidering, motivated by a mirror-photon exchange, refers 
to /n = 0. The allowed regions and bounds we will derive 
in the plane ti^^ - will therefore reflect this fact, and 
in particular will be shifted with respect to the standard 
case fp = fn which usually arises for many DM candi- 
dates. For analyses which relax the assumption fp — /„ 
see e.g. [HI Hi [Ml 133 ESHSS] . 

The long-range nature of the interaction occurs when 
iTi^ q. In this regime the differential cross section ac- 
quires an explicit dependence on the nuclear recoil energy 
and a Rutherford-like cross section emerges: 



da{v, Er) 
dE„ 



X 1 



2 F^Er) « E],^, 



(9) 



The E'^^ drop-off of the cross-section, which is not 
present in the point-like case, makes experiments with 
low energy thresholds (like DAMA and CoGeNT) to 
respond better to the interaction mechanism, as com- 
pared to experiments with relatively high energy thresh- 
old (CMDS and XENONIOO have stable thresholds of the 
order of /larger than 10 keV), and so may in principle im- 
prove the compatibility among those experiments. The 
recoil energy at which the interaction becomes effectively 
long-range depends on the mass of the target. In partic- 
ular, larger is the mass smaller is such transition energy. 



Therefore, for intermediate mass of the dark photon, a 
possible feature may arise, such that, at a given recoil 
energy, the interaction might be effectively long-range in 
experiments with large target-mass, while being in the 
point-like limit in low target-mass experiments. 

Considering typical a mass of targets tun ^ 100 GeV 
and nuclear recoil energy windows around few to tens of 
keV, the long range nature of the interaction manifests 
itself if the mass of the dark photon is lesser than 10 
MeV. 

More generally, the differential cross section of Eq. 
51) in terms of a normalized total cross section cr? 
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da{v,ER) _ ruN J_ ^2 p 



dER 



ml/{2mN) 



\ER + ml/i2mN) 



F'{Er), 
(10) 



where = 1 GeV ~ rup. We will use as free parameters 
in our analysis cr^^ and m^, an alternative choice being 
and the effective coupling of the DM with the nucleon 

As for the nuclear form factors, we adopt the standard 
form described by Helm in Ref. [M] : 

F(gr^) = 3^1^exp[-(gs)V2] (11) 
qrN 

where ji is the spherical Bessel function of the first kind 
with n = 1: 



sin(grAr) cos((jrAr) 



{qr 



N) 



(12) 



In Eqs. ( 11|12 ) rjy is the nuclear radius and s is es- 
timated trying to reproduce the more accurate results 
obtained from numerical evaluation of the Fourier trans- 
form relative to a Fermi distribution of scattering centers. 
A good agreement is obtained for s ~ (197MeV)~^ and 
rjv = ((l/(164MeV)Ai/3)2 _ 5s2)i/2^ remind that 

this expression of spin-independent form factor is derived 
assuming a Fermi distribution for the nuclear charge and 
that all the parameters used in this parameterization may 
be affected by sizable uncertainties. 

B. Rate of nuclear recoils 

The differential recoil rate of a detector can be defined 



as: 



dR f daiv,ER) 
dER -''^J ^iE^'"'''-^ 

where Nt is the total number of targets in the detec- 
tor {Na is the Avogadro's number) and dn^ is the local 
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FIG. 1: In the plane e vs. m^, iso-contours of constant rate 
(chosen as 1 cpd/kg) on a Na target, for a 10 GeV DM particle 
scattering and for various values for the energy threshold are 
shown. The galactic halo model is an isothermal sphere with 
Maxwell-Boltzmann (MB) velocity distribution with velocity 
dispersion vo ~ 220 km s"^ and local density po = 0.3 GeV 
cm"^. 



number density of DM particles with velocities in the el- 
emental volume (Pv around v. This last term can be ex- 
pressed as function of the DM velocity distribution (v) 
in the Earth frame, which is related to the DM velocity 
distribution in the galactic frame fciw) by a galilean ve- 
locity transformation as fE{v) — fci^ + i^Eit)) where 
vsit) is the time-dependent Earth (or detector) velocity 
with respect to the galactic frame. The prominent time- 
dependence (on the time-scale of an experiment) comes 
from the annual rotation of the Earth around the Sun, 
which is the origin of the annual modulation effect of the 
direct detection rate [T]. More specifically: 



VE{t) ^vg + vs + vmit) 



(14) 



The galactic rotational velocity of our local system vq 
and the Sun's proper motion vs are basically aligned and 
their absolute values are vq = vq ^ 220 ± 50 km s^^ and 
Vs = 12 km s~^, while the Earth rotational velocity V(^{t) 
has a size = 30 km s~^, period of 1 year and phase 
such that it is aligned to vq around June 2nd and it is 
inclined of an angle 7 ~ 60 degrees with respect to the 
galactic plane. More details can be found, for instance. 



in Ref. [65]. Summarizing: 



(15) 



where — ^^po/m^ is the local DM number density in 
the Galaxy and is determined by the local dark DM mat- 
ter density po and, in general, on a scaling factor ^ which 
accounts for the possibility that the specific DM candi- 
date under consideration does not represent the whole 
amount of DM. For all practical purposes, we can just 
assume ^ = 1 here. In Eq. ( 15 ) the velocity distribution 
function needs to be properly normalized. This can be 
achieved by requiring that in the galactic frame: 



d^foiv) = 1 



(16) 



where Vcsc denotes the escape velocity of DM particles 
in the Milky Way. For definiteness, we will adopt here 
Vesc = 650 km s^^. 

When considering the generalization of the differential 
cross-section to the case of long-range forces, the differ- 
ential rate of nuclear recoil can be cast in the form: 



dER 
where: 



2p^p 



G{Er) 



ml/{2mN) 



(17) 



F^iEn) (18) 



and: 



-^(^miii! ^) 



fEiv) 



d\ 



(19) 



with finin(£'i?) given in Eq. ([s]). The detection rate is 
function of time through the function X(t)i„in, t), as a con- 
sequence of the annual motion of the Earth around the 
Sun. The actual form of the function X{vnan,t) depends 
on the velocity distribution function of the DM particles 
in the halo. We will consider two cases: an isothermal 
sphere, whose velocity distribution function in the galac- 
tic frame fci^) is a Maxwell-Boltzmann function, and a 
triaxial model, with an anisotropic fci'^)- We will discuss 
in more details our choices in the next Section. 

Since the Earth velocity v^{t) has an explicit depen- 
dence on time due to the movement of the Earth around 
the Sun, and since this last velocity component is rela- 
tively small when compared to the main boost compo- 
nent represented by vq + vs, it is convenient to define a 
time-dependent parameter r]{t) as: 



r)E{t) ^T]Q+Ar] cos [27r(t - 0)/r] 



(20) 
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where ry© = {vg + vs)/vq and Arj = ti0COS7/vo, with 
A77 <^ ijQ, and where (p = 152.5 days (June 2nd) is the 
phase and r = 365 days is the period of the Earth mo- 
tion around the Sun. In Eq. ( 20 1 the time t is clearly 



Point-like {m^ = 1 GeV) 



expressed in days. By means of the approximation in Eq. 



( 20 1 we can define a convenient expansion of the recoil 
rate, which is suitable for velocity distribution functions 
which are not strongly anisotropic: 



dR 

dER 

dR 



it) 



dE 



R 



nE=ve 



d dR 



(21) 

Ajycos [27r(i - 0)/t] . 



Ve=Vq 



In order to properly reproduce the experimental recoil 
rate, we should account for effects associated with the 
detector response. We therefore need to include both 
the effect of partial recollection of the released energy 
(quenching) and the energy resolution of the detector. 
This can be done by the following energy-transformation 
and convolution: 





A 1\ /T A r^UT7 

JJAiVlA Litir 


LiO(si Mod 


MB Halo (t;o = 170 km/s) 


9.56 a 


1.90 a 


MB Halo (no = 220 km/s) 


9.57 a 


1.77 a 


MB Halo (t;o = 270 km/s) 


8.67 a 


1.46 a 


Triaxial Halo 


9.55 a 


1.82 a 



TABLE I: Statistical evidence of presence of modulation for 
the DAMA and CoGeNT data sets in the case of a point-like 
interaction (m^ = 1 GeV). This table refers to the analy- 
sis performed under the assumption of an isothermal sphere 
with a Maxwell-Boltzmann (MB) velocity distribution func- 
tion and for the triaxial halo model discussed in the text. 
For the MB case, the results for three different values of the 
dispersion velocity are reported. The analysis for the DAMA 
experiment refers to the annual modulation data, bounded by 
the total (unmodulated) rate I31|31[73], with the fraction of 
channeling varied in its allowed interval [HS]. The CoGeNT 
analysis considers the fit to the modulation amplitude [S], 
bounded by the total (unmodulated) rate with L-peaks sub- 
tracted. 



-|^(£;det) = / dE'IC{E,^,,E')J2§^ (eu = - 
dEdet J ^ dEji \ Qi 

(22) 

where the index i runs over the different nuclear species 
of the detector, E^ct is the detected energy and qi are 
the quenching factors of each of the nuclear species. The 
function /C(£'dct, E') takes into account the response and 
energy resolution of the detector, for which we assume, 
as is generally done, a gaussian behavior. 

As a final step, we need to average the recoil rate of Eq. 
(22 1 over the energy bins of the detector. We therefore 
define the unmodulated components of the rate 5ofe and 
the modulation amplitudes Smk for each energy bin k of 
width Ai?fc as: 



(m^ = 10 MeV) 



s„ 



Sok — 



1 



^ f dE 

AEk Jae.. dEdct 



(23) 



m=ne 



dE. 



d dR 



dct 



dm dE, 



dot 



At]. (24) 



Sok and Smk are the relevant quantities we will use in 
the analysis of the experimental data of DAMA and Co- 
GeNT. In the case of experiments which do not address 
the annual modulation effect, only Sok are relevant. 



III. DATASETS AND ANALYSIS TECHNIQUE 

Let us now move to the discussion of the techniques 
we use to analyze the various data sets. 





DAMA ChF 


CoG Mod 


MB Halo {vo = 170 km/s) 


9.17 a 


1.69 a 


MB Halo {vo = 220 km/s) 


8.77 a 


1.55 a 


MB Halo {vo = 270 km/s) 


7.74 a 


1.29 a 


Triaxial Halo 


9.01 a 


1.63 a 



TABLE II: The same as in Tab. |I] for the case of a mediator 
of mass = 10 MeV. 



For DAMA, CoGeNT and CRESST, we adopt the 
same technique of Ref. [BJ: we test the null hypoth- 
esis (absence of modulation for DAMA and CoGeNT, 
and absence of signal on top of estimated background 
for CRESST). From this test we obtain two pieces of in- 
formation: i) the level at which each data set allows to 
reject the null hypothesis (we will find a confidence level 
of about 8-9 cr for DAMA, 1-2 a for CoGeNT, and 4 
a for CRESST); ii) in the relevant DM parameter space 
(defined by the DM mass rriy^ and the DM-proton cross 
section c^^) we will determine the domains where the val- 
ues of the likelihood-function differ more than na from 
the null hypothesis (absence of modulation), and thus 
the corresponding evidence of the DM signal. We will 
use 71 = 7,8, n = 1, and n = 3,4 for DAMA, CoGeNT, 
and CRESST, respectively |6j. This choice (test of the 
null hypothesis) allows more proper comparison between 
the results arising from experimental data sets with dif- 
ferent statistical significances and, for the case of DAMA, 
allows to implement a requirement of a very high C.L. 
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(m^ = 30 MeV) 





DAMA ChF 


CoG Mod 


MB Halo {vo = 170 km/s) 


9.50 a 


1.84 a 


MB Halo (vo = 220 km/s) 


9.40 a 


1.69 a 


MB Halo {vo = 270 km/s) 


8.39 a 


1.40 cr 


Triaxial Halo 


9.52 (7 


1.77 a 



TABLE HI: The same as in Tab. ID for the case of a mediator 
of mass = 30 MeV. 



Long range (m^ = MeV) 





DAMA ChF 


CoG Mod 


MB Halo {vo = 170 km/s) 


8.76 a 


1.61 CT 


MB Halo {vo = 220 km/s) 


8.23 a 


1.48 a 


MB Halo {vo = 270 km/s) 


7.27 (7 


1.23 a 


Triaxial Halo 


8.64 a 


1.55 a 



TABLE IV: The same as in Tab. ID for the case of a long- 
range interaction (m^ — 0). 

Our statistical estimator is a likelihood function, de- 
fined as £ = Yli I where i stands for the i-th energy 
bin in DAMA and CoGeNT, and for the i-th detector in 
CRESST. Ci is the likehhood of detecting the number of 
observed events given the expected background and DM 
signal. Ci are taken to be gaussian for DAMA and Co- 
GeNT and poissonian for CRESST, since in this case the 
number of events in each sub-detector is low. Defining 
Cbg as the likelihood of absence of signal (i.e., without 
the DM contribution), the function y = — 21n£f,g/£ is 
assumed to be distributed as a x^~variable with one de- 
gree of freedom, for each value of the DM mass (note 
that in the DAMA and CoGeNT cases, y simply reduces 
to y ^ xlg — X^)- From the y function we extract, for 
each value of the DM mass, the interval on tr^^ where 
the null hypothesis (absence of modulation, i.e. a^^ = 0) 
can be excluded at the chosen level of confidence: 7a 
(outer region) or Scr (inner region) for DAMA and la 
for CoGeNT). From this, regions in the - cr^^ plane 
arise. 

Constraints from null experiments are derived by 
constructing again a similar likelihood-function A — 
—2\nC/£hg, where £ is the likehhood of detecting the 
number of observed events (2 and 3 for CDMS and 
XENONIOO, respectively) over the whole energy range 
of the experiment given the expected background and 
the DM signal, while in the DM signal is not in- 
cluded. Both likelihoods are taken to be Poissonian and 
A is assumed to follow a x^-distribution. Bounds are 
conservatively shown at 5-a C.L. 



A. DAMA 

DAMA, located at the Laboratori Nazionali del Gran 
Sasso, is an observatory for rare processes based on the 
developments and use of highly radiopure scintillators. 
The former DAMA/Nal and current DAM A/Libra ex- 
periments, made of radiopure Nal(Tl) crystals, have the 
main aim of investigating the presence of DM particles in 
the galactic halo by looking at their annual modulation 
signature. 

The signal in DAMA is the energy deposited in scintil- 
lation light. On the other hand, the scattered nucleus is 
losing energy both electromagnetically and through nu- 
clear interactions: this effect is taken into account by the 
quenching factors q which convert the total nuclear recoil 
energy to the energy seen by the detector i^dot = qE^- 
For NaT crystals we take, QNa = 0.3 and qi = 0.09 [57] . 
Notice that the uncertainty on the actual values of the 
quenching factors in NaT 'SS' can have a visible impact 
on the reconstructed DM properties [6]. 

However, it has been appreciated [HTJ EH] that nuclei 
recoiling along the characteristic planes of crystals can 
travel large distances without colliding with other nuclei, 
and essentially deposit all their energy electromagneti- 
cally (which corresponds to g = 1 and thus Edet = Eb,). 
This process is known as channeling and its relevance in 
DM direct detection experiments (and in particular for 
Nal crystals) is currently under scrutiny [69^171] . Since 
the actual amount of channeling in a detector like DAMA 
is not currently know, we take the following three-fold 
approach: i) we show the effect induced by a sizable chan- 
neling effect, at the level estimated by the DAMA Collab- 
oration in Ref. |69j . by employing an energy-dependent 
channeling fraction fch as reported in Fig. 4 of [69 ; ii) we 
consider the case of a negligible channeling effect; iii) due 
to these uncertainties in the knowledge of the amount 
of channeling, and considering that the actual value of 
fch is likely to lie between the two previous cases, we 
smoothly vary fch between them and marginalize over 
it. We consider this last approach as the most general 
and we will adopt it for most of our analyses. Notice 
that a small fraction of channeling is exactly what would 
easily reproduce a clear agreement between DAMA and 
CoGeNT results also in the standard case of a point- 
like interactions (see for instance the results in Ref. [5]). 
Therefore to allow for a variation of fch appears to be a 
useful approach, when comparing results from different 
experiments. 

We consider the whole set of DAMA/Nal [3| and 
DAM A/LIBRA 4J data, which correspond to a cumu- 
lative exposure of 1.17 tonxyr. We analyze the mod- 
ulation amplitudes S'^j^ reported in Fig. 6 in Ref. [4] 
by using our statistical technique discussed above. The 
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modulation amplitudes of Ref. [1] can be considered as 
a data-reduction of the time- and energy- dependent 
data in 8 energy bins. The actual values of of Ref. 
[3] are valid under the assumption that annual modu- 
lation occurs with phase and period fixed at day 152.5 
and 365 days, respectively. We do not use directly the 
time-series of the data (which would be a better option 
for our statistical technique of studying the null hypoth- 
esis) since these are available only in three energy bins 
in Rcfs. [3l Hj, while instead the S'^J^ are provided in 
8 energy bins [1], therefore supplying more information 
especially for light DM. We checked that the results ob- 
tained by using the S'^^ is in full agreement with the 
results of Refs. [6], where the same type of statistical 
analysis has been employed directly on the DAMA data. 

In addition, we use the information on the total rate 
as a constraint, by requiring that the DM contribution 
^0 does not exceed the corresponding experimental value 
Sq^^ in the 2-4 keV energy range, measured by DAMA 
[73] . Summarizing, for the DAMA datasets our approach 
requires to determine: 



y = -21n/: = x^(e, "1,0,771^) 

E jSmk — '5'toI') , {Sq - 
k=l ^ 



(25) 
0(^0-^0'^''), 



where and a are the experimental errors on 5"°^^ and 



S'q^^, respectively, and Q denotes the Heaviside function. 



The last term in Eq. ( 25 1 implements the upper bound on 



the unmodulated component of the rate 5*0, by penalizing 
the likelihood when 5*0 exceeds S"^^ . For the detector 
energy-resolution we use a Gaussian function of width 
a,^^{E) = £;(0.448//E-h 0.0091) I72j. 

B. CoGeNT 

The CoGeNT experiment is made by Ge detectors with 
very low threshold (about 0.4 keVee, where keVee de- 
notes keV electron-equivalent energy). Thanks to this 
property, CoGeNT has the capability of being very sen- 
sitive to DM candidates with rriy. < 10 GeV, although 
large background contamination may be present at these 
low-energies. 

In 2010, the CoGeNT Collaboration reported the de- 
tection of an excess not identifiable with a known back- 
ground, and potentially compatible with a DM interpre- 
tation 7 . More recently, the temporal evolution of the 
measured rate in different energy bins for data taken be- 
tween December 4, 2009 and March 6, 2011 (442 live- 
days) has been presented [5] . At low energies, CoGeNT 
data favour the presence of an annual modulation which 
can be fitted by a WIMP signal [S]. 



One could perform an analysis similar to the one real- 
ized for DAMA by deriving the modulation amplitudes 
S^^l from the data presented by the COGeNT Collab- 
oration in Fig. 4 of Ref. [8]. To this aim we fix the 
period of modulation at 365 days and the phase at June 
2nd (day 152.5), as for the DAMA dataset. We obtain 
S^^l = 0.12 ± 0.08 counts/day in the 0.5-0.9 keV energy 
bin and S"^^ = 0.26 ± 0.17 counts/day in the 0.9-3.0 
keV energy bin. However, as mentioned above, for our 
statistical technique (test of the null hypothesis) acting 
directly on the timeseries of the data is more appropriate, 
and so for CoGeNT we follow this path. In addition to 
the modulation amplitude, we treat the total rate mea- 
sure by CoGeNT as a constraint and, similar to the case 
of the analysis of the DAMA data, we define: 



y = -2\nC = x^(e, 777,0,771;^) 



(26) 



^ — -r^^^^ — ^k — 

k=l fc=l 
31 

E 



(5oj - SqY) 



2 ^{Soj - S"^^) 



where S^k = 1/Atfc /aj^. '5'™^ cos [27r(t - (?!))/r], with 
Aifc being the temporal bin of experimental data, and 
SZl = RZl- < RZ" > with RZl being the total rate 
shown in Fig. 4 of Ref. [S] and < > being the 

rate R'^^ averaged over a cycle (1 year). We compute 
the total rate in the 0.9 — 3.0 keV energy-bin by simply 
subtracting the rate in 0.5 — 0.9 keV bin to the rate in 
the 0.5 — 3.0 keV bin, and with a Gaussian propagation 
of errors. Sq^^ and aj denote the counts and their corre- 
sponding errors given in Ref. |B] (31 bins in the energy 
interval 0.4 — 2 keVee), with the L-shell peaks removed, 
but without any further background removal. 



Notice that with the procedure of Eq. (261 we do not 



require that the whole total (unmodulated) spectrum 
of CoGeNT is due to DM scattering: we just require 
that the the spectrum is not exceeded by our theoretical 
model. The unmodulated spectrum acts therefore as 
a bound, leaving room in it for an unknown background 
component. 

The total fiducial mass of the CoGeNT experiment is 
330 g, the energy resolution is described by a Gaussian 
function with the form of the width a^cs taken from |74| , 
and the quenching factor follows from the relation E = 
0.2 E^^^ below 10 keV [75j. 
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C. CRESST 

The CRESST experimental setup [ST] at the Gran 
Sasso Laboratories includes 300 g of scintillating CaW04 
target crystals. The particle interaction is detected 
through phonons by the phase transition thermome- 
ter and through scintillation light by a separate cryo- 
genic detector. Results from 730 kg days of data-taking 
have been recently presented by the CRESST collabo- 
ration [SH]. Sixty-seven events are found in the WIMP 
acceptance region, and background contributions from 
leakage of e/7-events, neutrons, a-particles, and recoil- 
ing nuclei in a-decays are not sufficient to explain all 
the observed events. A likelihood-ratio test rejects the 
background-only hypothesis at a significance larger than 
4-a [55]. 

To perform our analysis, we compute DM signal-events 
in each of the 8 CRESST detector modules. We consider 
acceptance regions and number of observed events, as 
provided in Table 1 of Ref. [58 , and background events 
are derived according to estimates in Sec. 4 of Ref. [58] . 
Performing a likelihood-ratio test, we obtain an evidence 
for the best -fit DM signal over a background-only sce- 
nario at 4.1-(T C.L., thus in good agreement with the 
result of the Collaboration. 

However, we do not have at disposal all the information 
needed to perform a full analysis, so our derived contours 
and allowed regions must then be considered as indica- 
tive. In particular we use the published exposure of 730 
kg X days assuming an equal contribution from each mod- 
ule (which thus has an exposure of 730/8 kg x days) and 
consider a constant efficiency. Moreover, in order to be 
able to properly discriminate among DM models induc- 
ing different recoil spectra, we would need the energies of 
events in each detector, rather than the total number in 
the whole acceptance region. This is the reason why our 
allowed regions for the standard point-like case do not 
overlap with the regions presented in Ref. [35] for large 
dark matter masses. On the other hand, at masses be- 
low 50 GeV (which is the main focus of this paper), the 
agreement becomes very good. 

The likelihood function y is used as described above to 
determine the allowed regions in the parameter space of 
our model. 

D. CDMS 

The CDMS experiment, located at the Soudan Un- 
derground Laboratory in Minnesota (like CoGeNT), has 
operated cryogenic semi-conductor detectors to simulta- 
neously measure phonons and charge in order to reject 
most of the dominant radioactive background and dis- 



entangle a DM signal. Only the Ge detector were com- 
pletely functioning and fully exploited to set constraints 
on DM properties. 

In this paper, we consider the 'standard' 2009 CDMS- 
II results They are based on Ge data taken between 
July 2007 and September 2008, applying conservative nu- 
clear recoil selection cuts and assuming a 10 keV energy 
threshold. The total exposure is 612 kgx days and ef- 
ficiency has been taken from Fig. 5 in Ref. [75] (black 
curve) with quenching factor g ~ 1.^ Two signal events 
were seen in the 10-100 keV energy window [SS] against 
an expected background of 0.9 ±0.2 (which are the num- 
bers we use to derive constraints). 



E. XENONIOO 

The XENONIOO experiment searches for DM scatter- 
ings on a target of purified liquid Xe by measuring scin- 
tillation and ionization signals. We consider results pre- 
sented in Ref. [60], corresponding to an exposure of 100.9 
days in a 48 kg fiducial volume. After the implementa- 
tion of all the cuts three events have been reported in 
the DM signal region with an expected background of 
1.8 ±0.6 events. 

The XENON collaboration selected the energy win- 
dow for the WIMP search region to be between 4-30 
photoelectrons (in terms of prompt scintillation light in 
the liquid), corresponding to recoil energy of 8.4 - 44.6 
kcVnr (based on the C^s parametrization in Ref. j60|). 
However, recoils at lower energy can contribute as well, 
especially close to threshold, due to the Poissonian tail. 
Both the statistics and the quenching of few keVnr nu- 
clear recoils in liquid Xenon are not completely under- 
stood (for recent discussions, see e.g., Ref. [52] and ref- 
erence therein). For definiteness, we consider a Poisso- 
nian distribution of photoelectrons, and include a single- 
photoelectron resolution of 0.5. The £off function is a 
very crucial ingredient at this low level of photoelectrons 
and for light DM. We adopt two different approaches, in 
order to bracket a possible (but definitely not exhaustive) 
uncertainty on the derived bounds from XENONIOO: i) 
we adopt the nominal central value of £off shown in Fig.l 
of Ref. [60] , which heavily relies on linear extrapolation 
below 3 keVnr; ii) more conservatively, we increase the 
photomultiplier threshold to 8 photoelectrons, in order 
to determine a situation which is nearly independent on 



For the cases of interests (low DM masses), a similar analysis 
can be performed by exploiting combined data from CDMS and 
EDELWEISS experiments (see Fig. 1 in Ref. 77 ), obtaining 
basically the same results. 
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the knowledge of £eff below 3 keVnr. The value of 8 pho- 
toelectrons has been chosen to this purpose (namely, it 
is the lowest value satisfying such requirement). Notice 
that these two recipes do not exhaust the possibilities of 
alternative assumptions that can be done to determine 
the XENON 100 response to light DM. For more discus- 
sion and additional considerations, see e.g. Refs. [STll62j . 
Due to the large uncertainties inherent in the derivation 
of bounds from XENONIOO for hght DM, it appears to 
be still preliminary to assume those bounds as strictly 
firm. We nevertheless show them, but without enforcing 
them in our discussion. In any case, conservatively, we 
consider the 8 PE bound as more appropriate since it is 
less dependent on the Lcs extrapolation. 

Finally, to compute the expected signal we follow Eq. 
(13-16) in Ref. In both the cases of CDMS and 

XENON, we derive upper bounds as discussed at the 
beginning of Section |III[ 

IV. RESULTS 

Let us now move to the discussion of our results in 
terms of long-range interactions of DM on protons. We 
show in this Section the effect of moving from the point- 
like case (parameterized, for convenience, in terms of a 
massive mediator (/), for which we use a reference value 
for its mass of 1 GcV, larger values reproducing the same 
results) to the very long range case (which refers to the 
extreme case of = 0) , going through the intermediate 
case of relatively light <j). To this aim, we sit where the 
sensitivity to is largest, which is for values ranging 
from several MeV's to few tens of MeV's. This fact can 
be seen, in an illustrative case, in Fig. [T] where we show, 
in the plane e vs. m^, the iso-contours of constant rate 
(chosen as 1 cpd/kg) on a Na target, for a 10 GeV DM- 
particle scattering and for various values for the energy 
threshold of a Na detector (for definiteness, we set here 
Z' — 1 and Q!dark = c^sm)- As discussed in Section [lH 
for values of the mediator mass smaller that a few MeV 
(the actual value depending on the recoiling nucleus and 
detected energies) the differential cross section and the 
corresponding experimental rate become proportional to 
e'^/Ej^, typical of a very long-range interaction, and the 
dependence of the iso-contours of constant rate becomes 
effectively insensitive to m^. Fig. [Tjshows that this long- 
range limit is reached, for typical nuclei and recoil ener- 
gies in the range of interest of actual experiments, for 
7710 < 1 MeV. On the contrary, for values of larger 
that several tens of MeV, the cross section reaches the 
point -like limit: in this case the detection rate is pro- 
portional to e^/m^ and the iso-contours of constant rate 
follow the linear behaviour (in log-log scale) shown in the 



figure. The transition regime is then obtained for for 
values ranging from few MeV to few tens of MeV. For 
definiteness, in our analyses we will adopt the values of 
TO^ = 10 MeV and 30 MeV. 

Since the response of direct detection experiments is 
quite sensitive to the DM distribution in the galactic halo 
[5], especially in velocity space, we perform our analysis 
by considering two kinds of variations in this respect. 

First, we analyze the direct detection datasets by us- 
ing a standard isothermal model, which basically im- 
plies a truncated Maxwell-Boltzmann velocity distribu- 
tion function, but we take into account uncertainties on 
the velocity dispersion ?;o, as discussed in Ref. [5]. We 
recall that in the case of an isothermal model, the DM ve- 
locity dispersion is directly linked to the asymptotic value 
of the rotational velocity supported by the DM halo: un- 
certainties in the velocity dispersion vq are therefore rep- 
resentative of the uncertainties in the local rotational ve- 
locity 0. Following Ref. |5] we will use the three values 
Vq — 170, 220, 270 km s~^, which bracket the uncertainty 
in the local rotational velocity. Let us notice that also 
the value of the local DM density pq is correlated to the 
adopted value of vq, as discussed e.g. in Ref. [S]. This 
case corresponds to the model denoted as AO in Ref. [S], 
and we adopt the case of minimal halo, which implies 
lower values of local DM density (since a fraction of the 
galactic potential is supported by the disk/bulge). In 
turn, this implies the adoption of po = 0.18,0.30,0.45 
GeV cm"3 for vq = 170,220,270 km s"\ respectively 

Second, we consider a different form for the galac- 
tic halo, and in order to somehow emphasize the dif- 
ference with the isothermal sphere, we adopt a triaxial 
halo model 79J, with anisotropic velocity dispersions. We 
adopt the model denoted by D2 in Ref. [5], which corre- 
sponds to the case when the Earth is located on the major 
axis of the potential ellipsoids, with vq = 220 km s^^ and 
maximal halo, which implies po = 0.84 GcV cm^'^ 0. 

A. Point— like scenarios 

We start with the point -like scenario. Fig. [2] shows 
the scattering cross sections on proton cr^^ as a func- 
tion of the dark matter mass m^. The galactic halo 
has been assumed in the form of an isothermal sphere 
with velocity dispersion vq — 220 km s~^ and local den- 
sity po = 0.3 GcV cni^-^. In the left panel wc show the 
allowed regions compatible with the annual modulation 
effects in DAMA, the CoGeNT excess, as well as the re- 
gion which turns out to be compatible with the CRESST 
excess, when interpreted as a DM signal. Specifically, 
the solid green contours A, denote the regions compat- 
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FIG. 2: Point-like scattering cross sections on proton, as a function of the dark matter mass. The galactic halo has been 
assumed in the form of an isothermal sphere with velocity dispersion vo ~ 220 km s~^ and local density po = 0.3 GeV cm~^. 
Left panel: The solid green contours A, denote the regions compatible with the DAMA annual modulation effect |3l|4], in 
absence of channeling. The solid red contours B, refer to the regions compatible with the DAMA annual modulation effect, 
when the channeling effect is considered at its maximal value. The dotted blue contour refers to the region derived from the 
CoGeNT annual modulation effect 8 , when the bound from the unmodulated CoGeNT data is included. The dashed brown 
contours denote the regions compatible with the CRESST excess |58) . For all the data sets, the contours refer to regions where 
the absence of modulation can be excluded with a C.L. of: 7a (outer region), 8a (inner region) for DAMA, la for CoGeNT and 
3a (outer region), ia (inner region) for CRESST. Right panel: The same as in the left panel, with the following difference: 
the solid orange contour refers to the DAMA annual modulation data, when the fraction of channeling is varied in its allowed 
interval [69]. Again, the contour refers to the region where absence of modulation can be excluded with a C.L. 7a, Sa. 



ible with the DAMA annual modulation effect [3j |4] , in 
absence of channeling. The solid red contours B, refer to 
the regions compatible with the DAMA annual modula- 
tion effect, when the channeling effect is considered at its 
maximal value. ^ The dotted blue contour refers to the 
region derived from the CoGeNT annual modulation ef- 
fect ^ , when the bound from the unmodulated CoGeNT 
data is included. The dashed brown contours denote the 
regions compatible with the CRESST excess (SSj. For 
all the data sets, the contours refer to regions where the 
absence of modulation can be excluded with a C.L. of: 
7(7 (outer region), 8a (inner region) for DAMA, la for 
CoGeNT and 3cr (outer region), Aa (inner region) for 



^ DAMA is a multi-target detector and allowed regions at large 
DM mass correspond to scattering on I, while at small DM mass 
regions correspond to scattering on Na. Note that the region at 
DM mass around 10 GeV is given by scattering on Na targets 
in the no-channeling case, and by scattering on I targets in the 
maximal channeling case. 



CRESST. The right panel, instead, shows a further anal- 
yses for the DAMA data: the solid orange contours refer 
to the results obtained by varying the channeling fraction 
fch in its allowed range, as discussed in Section pll A[ We 
can therefore see the extent of the DAMA allowed region 
when /ch is marginalized over. We will adopt this pro- 
cedure of treating the channeling effect in DAMA in the 
remainder of the paper. 

From Fig. [2] we notice that, in the case of a point- 
like cross section, the DAMA and CoGeNT regions both 
point toward a DM with a mass in the 10 GeV ballpark 
(more specifically, from about 5 up to about 20 GeV) and 
cross sections from a few x 10"^-'^ cm^ to 10"'^^ cm^ (with 
our choice of galactic halo model and for our situation 
of scattering on protons only), a result which confirms 
similar analyses [SlIMT^. The DAMA and CoGeNT re- 
gions are largely overlapped. This is especially clear in 
the right panel of Fig. [2] where we marginalize over the 
channeling fraction: even a small amount of channeling 
is enough to make DAMA and CoGeNT regions perfectly 
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FIG. 3: Point-like scattering cross sections on proton, as a function of the dark matter mass. The galactic halo has been 
assumed in the form of an isothermal sphere with velocity dispersion i;o = 170 km s~^ and local density po = 0.18 GeV cm~^ 
(left panel); vo ~ 270 km s~^ and local density po = 0.45 GeV cm~^ (right panel). Notations are the same as in the left panel 
of Fig. [2] 



overlap. 

It is interesting to note that, when channeling is in- 
troduced, there is an allowed region for DM mass of 4-5 
GeV, which corresponds to the signal given by the scat- 
tering of DM particles on the fraction of Na targets un- 
dergoing channeling. This region is very seldom consid- 
ered in analyses of DAMA data, and the reason resides in 
the statistical technique adopted. Indeed, the statistical 
evidence associated to this region is very high when com- 
pared to the absence of signal (up to 8-a), but its is 
substantially larger than the best-fit DM model (which 
is for DM mass around 10 GeV). Therefore, when regions 
are drawn assuming that a DM signal is present and in- 
cluding models that falls within a certain C.L. from the 
best-fit case (as often done in the literature, except in the 
analysis performed by the DAMA Collaboration itself, 
see e.g. Refs. [Ml 150] and Refs. [5il6J) such region disap- 
pears (unless allowing extremely large C.L. with respect 
to the minimum x^)- This is also the case in our Fig. 13, 
where we show "preferred" regions, i.e., domains which 
include models having a within a given C.L. with re- 
spect to the best-fit x^- An allowed region for DM mass 
of 4-5 GeV is particularly interesting, especially in light 
of the fact that it is compatible with both DAMA and 
CoGeNT data, and can easily satisfy constraints from 
other experiments, such CDMS and XENONIOO (inde- 
pendently on the method employed to derive the latter) . 



The statistical significance for the presence of a signal 
due to DM scattering is reported in Table |lj which con- 
firms that for DAMA the interpretation in terms of a DM 
signal is highly favored {9.6a for the case under consid- 
eration, i.e. isothermal sphere with central value of vq). 
In the case of CoGeNT the effect is at the level of l.Str, 
a fact that simply reflects the current lower statistics of 
the CoGeNT data sample. 

Let us stress that we are not considering here a number 
of sources of uncertainties of experimental origin which 
can be relevant in the analysis, like for instance the un- 
certainties on the quenching factors in Na, I and Ge [S]. 
Those effects would somehow enlarge both the DAMA 
and CoGeNT (and CRESST, too) regions, allowing a 
wider range of DM masses. 

Fig. [2] also shows, as brown contours, the preferred re- 
gion derived from the CRESST excess. We remind that 
some assumptions on unknown features of the data were 
needed. In particular, without energy spectra of recoil 
events for each module, the discrimination between a 
large-mass WIMP (nearly flat spectrum) and a small- 
mass WIMP (exponential spectrum) becomes difficult. 
The total spectrum shown in Fig. 11 of Ref. [SS] points 
toward the latter rather than the former case. There- 
fore, at small masses, we do not expect CRESST regions 
to be drastically modified by an analysis with the full 
set of information. Indeed, our preferred region perfectly 
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FIG. 4: Scattering cross sections on proton, as a function of the dark matter mass, for a long-range mediator of mass = 10 
MeV (left panel) and 30 MeV (right panel). The galactic halo has been assumed in the form of an isothermal sphere with 
velocity dispersion vq = 220 km s~^ and local density po — 0.3 GeV cm~^. Notations are the same as in the left panel of Fig. 

m 



overlaps with the one found in [58] for DM mass below 50 
GeV. The picture is different at large masses, where our 
CRESST regions are affected by our imperfect knowledge 
on relevant detection information. 

In the current situation, the analysis shows that the 
contours perfectly overlap with the DAMA and CoGeNT 
regions, and it is very intriguing that all three "positive" 
experimental results point to the same sector of DM pa- 
rameter space. 

The effect induced by the variation in the DM disper- 
sion velocity is shown in the two panels of Fig. |3] In 
the case of vq = 170 km s^^, regions are not significantly 
modified as compared to the case of vq = 220 km s~^ (ex- 
cept for the overall normalization, a fact which reflects 
the different values of local DM density in the two cases) . 
Table |T] shows that also the statistical significance is simi- 
lar to the Vq = 220 km s~^ case. Larger dispersion veloc- 
ities, instead, are more favorable for lighter DM masses. 
The right panel of Fig. [3] shows in fact that the allowed 
regions are shifted toward lower masses for vq = 270 km 
s"^. Regions still overlap, although they shrink to more 
defined ranges, both in DM mass and cross section. In 
this case, the DM mass cannot extend much above 10 
GeV. From Table|T]we see that the statistical significance 
is reduced by about la for DAMA. Clearly this level of 
reduction is not sizable enough to allow exclusion of the 
vo = 270 km case. 



B. Long— range forces scenarios 

We now proceed to discuss the situation when the 
mediator mass falls in the range of transition from 
the point-like to the non point-like case. Let us first 
summarize the effect of long-range interactions in di- 
rect detection experiments. As described in Section |ITj 
da/dEn cx {Eu^ + m'^/{2mN))~^ , which teUs us that 
larger is the target mass mjv smaller is the recoil en- 
ergy corresponding to the transition from the point-like 
to the non point-like case, and that, in the long-range 
limit, events at low recoil energy are enhanced with re- 
spect to the ones occurring at large recoil energy. Those 
simple arguments can be applied to understand the be- 
haviour of the unmodulated signal. For the modulated 
signal, on the other hand, the picture is less straightfor- 
ward. Indeed, if the minimal velocity providing a nu- 
clear recoil can be significantly smaller than then the 
(cosine-like) modulated signal is suppressed. It can hap- 
pen at low recoil energy (see Eq. [s]), and in this case the 
long-range interaction would not enhance the modulated 
signal. 

A first set of results are shown in Fig. [4j where 
we report the results for the two reference values of 
— 10, 30 MeV. The galactic halo is in the form of an 
isothermal sphere with local dispersion velocity vq = 220 
km s^^ and local density po = 0.3 GeV cm""^. The panels 
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FIG. 5: The same as in Fig. [4] for an isothermal sphere with velocity dispersion vo = 170 km s~^ and local density po — 0.18 
GeV cm-^ 
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FIG. 6: The same as in Fig. [4] for an isothermal sphere with velocity dispersion vo = 270 km s ^ and local density po ~ 0.45 
GeV cm-^ 



on Fig. |4]can be directly compared to the right panel of 
Fig. [2j We notice that moving from the point-like case to 
the case of an intermediate-mass mediator, the allowed 
regions are sizeably modified (especially for DAMA). In 
the case of DAMA, the DM low-mass region becomes in- 
creasingly preferred over the 100 GeV DM mass region, 



which basically disappears for — 10 MeV. The statis- 
tical significance of this cases is similar to the point-like 
case, as can be seen in Tables |ll] and |III| The same oc- 
curs in the case of CoGeNT. Following the arguments 
explained above, we find that in the case of scattering 
on I in DAMA, the unmodulated signal is enhanced at 
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FIG. 7; The same as in Fig. |4] except that the interaction 
is point-like, and for a triaxial galactic halo with the Earth 
located on the major axis [79], with velocity dispersion vq = 
220 km s"^ and local density po = 0.84 GeV cm"^ [5]. 



low energy by the long-range interaction while the cor- 
responding modulated signal is not, and so the region at 
mass above 50 GeV (and most of the channeled region 
at 10 GeV) is ruled out. Note also that the enhance- 
ment of the signal on Na at very low recoil energies cuts 
out the region with large cross-section and very small 
masses. In the CRESST case, we consider a signal in O, 
Ga and W, with the target-mass of the formers being 
significantly smaller than the latter. This means that, in 
the W case, the long-range interaction nature manifests 
itself at larger mediator mass (or, equivalently, at larger 
recoil energy). Therefore, while in the point-like case 
the region at large DM masses is mostly given by signal 
on W, the contributions of O and Ca becomes increas- 
ingly important as the mediator mass decreases (strictly 
speaking, when the mass of the mediator decreases, also 
the cross section needed to reproduce a given signal de- 
creases, until the long-range limit is reached, see Fig. 
1, but this process stops earlier in the W case). In the 
GoGeNT case, the region remains basically unchanged, 
partly because we have at disposal only two energy bins, 
a fact that makes hard any disentanglement of energy- 
dependent effects. 

In order to investigate the dependence of the results 
on some astrophysical assumptions, we report in Fig. [s] 
and [g] the analysis for intermediate-mass mediator (same 
reference values: = 10, 30 MeV) when the dispersion 



velocity of the MB distribution is changed: Fig. [5] refers 
to vq — 170 km s~^ (which, for consistency in the halo 
modeling required a local density po = 0.18 GeV cm~^ 
[5]), while Fig, |6] shows the case of vq — 270 km s~^ 
(po = 0.45 GeV cm~'^ [S]). As expected from kinematics 
of the DM scattering process, lower average velocities of 
DM (smaller values of vq) require larger values of the DM 
mass in order to reproduce the same effect in a detector: 
this implies an extension of the allowed regions toward 
larger DM masses for va = 170 km s^^. This is man- 
ifest in Fig. [5j The opposite is true when vq becomes 
large, as is clear from Fig. [6) Tables [lT| and |III| show that 
all these models with intermediate-mass mediators are 
viable at the same level as the point-like scenario, with 
no significant variation of the statistical significance, ex- 
cept again for the case vq = 270 km which gets worse 
when the mass of the light mediator becomes small. This 
is confirmed in Table |III[ where the case = is re- 
ported. We can conclude that large dispersion velocities 
and light mediators are slightly disfavored by the DAMA 
data (still not at the level of considering them as ex- 
cluded), while in the case of dispersion velocities in the 
range from vq = 170 km to about vq — 220 - 250 km 
s^^ the statistical agreement is basically stable over vari- 
ation of and Vq. In the case of CoGeNT, we do not 
observe significant variation of the statistical significance 
when we change m^: mild preference toward smaller val- 
ues of Vq are present also with intermediate-mass medi- 
ators, like in the case of point-like interactions, but not 
statistically significant. 

To further discuss the dependence on astrophysics, we 
extend the analysis to a different velocity distribution 
function, by adopting an anisotropic halo model, instead 
of the isotropic MB velocity distribution. For this, we 
adopt a triaxial halo model where the Earth is located 
on the major axis and the velocity anisotropy is tangen- 
tial [70]. The anisotropic velocity dispersion is taken as 
in model D2 of Ref. [5]. The results for this class of 
halo models are shown in Figs. [7] and [8] for the case of 
a point-like cross section (Fig. [t] ) and in the case of 
mediator masses = 10,30 MeV (Fig. [s]). This model 
presents a relatively large degree of anisotropy in the ve- 
locity dispersion, and this sizably modify the recoil rate 
[5] . When compared to the isotropic MB case, the triaxial 
halo model tends to somehow enlarge the range of allowed 
DM masses, especially in the case of intermediate-mass 
mediator. From the Tables, we notice that the triaxial 
halo models are all viable. 

The case of very long-range forces is shown in Fig. [9] 
and [To] for the isothermal sphere and in Fig. [TT] for the 
triaxial halo model, where we set the extreme case of 
= 0, and the statistical significances are reported in 
Table IV We notice that in the case of = 0, the al- 
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FIG. 8: The same as in Fig. [4] for a triaxial galactic halo with the Earth located on the major axis [79], with velocity dispersion 
vo = 220 km s~^ and local density po = 0.84 GeV cm~^ |5]. 
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for a dispersion velocity vq = 220 km s~^ and local density 
po = 0.3 GeV cm-^ 



lowed regions moves toward lighter DM, even sizably like 
in the case of an isothermal sphere with large dispersion 
velocities. However this last case is disfavored by the 



statistical analysis shown in Table |IV[ while the analysis 
of the CoGeNT data do not exhibit significant variation 
with the change of halo model, in the case of DAMA a 
very long-range interaction together with large disper- 
sion velocities is significantly worse (by about 2(t) than 
the other cases under analysis. The CoGeNT, CRESST, 
and DAMA data can have a common DM interpretation 
and the compatibility increases for large-intermediate 
and small-intermediate Vq (see Figs. [2 11), leading 
to similar conclusions as the ones derived above for each 
single datasets. 

Constraints that can be derived from the null experi- 



ments (CDMS and XENON 100) are shown in Fig. [12 
The different lines refer to the various galactic halo mod- 
els discusses in our analyses: broken lines refer to the 
isothermal sphere with vq — 170 km s~^ and pq — 0.18 
GeV cm^^ (short-dashed line), ?;o = 220 km s~^ po = 0.3 
GeV cm~^ (medium-dashed line), fo = 270 km and 
Po — 0.45 GeV cm"'^ (long-dashed line). Solid lines refer 
to the triaxial halo model [51 [75] . For the XENON de- 
tector, all the constraints refer to a threshold of 8 photo- 
electron. In the first panel, the blue dashed line stands, 
instead, for a threshold of 4 photoelectron and for an 
isothermal sphere with vq = 220 km s^^ and po = 0.3 
GeV cm~^. We can notice the extent of variation of 
the constraints when the galactic halo model and/or the 
mechanism of interaction is varied. As expected, when 
comparing long-range interactions with the point-like 
case, the relative impact of the bounds from XENON 
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FIG. 11: The same as in Fig. [4] except that the interaction 
is long-range (m^ — MeV), for the triaxial halo model. 



and CDMS at small DM mass is enhanced with respect 
to large mass, as can be seen by comparing the first panel 
to the other three. However, this effect is less pronounced 
than in experiments with lower threshold (such as DAMA 



and CoGeNT). 



C. Constraints on the mass of the light mediator 
and on the kinetic mixing 

We finally move to derive bounds on the light-mediator 
mass 7710 and on the kinetic mixing parameter e, under 
the hypothesis that the mechanism of long-range forces is 
compatible with the annual modulation results in direct 
detection. 

As a first analysis, we discuss the bounds on the me- 
diator mass by showing the allowed regions in the plane 
777,0 vs. the dark matter mass 777^, obtained by adopting a 
maximum likelihood method. For definiteness, we adopt 
the DAMA data set, since this shows a clear modulation 
effect with a very large C.L. The results are reported in 
Fig. [13] The case of an isothermal sphere with vq = 220 
km s^^ is shown in the left panel, while the triaxial halo 
case is reported in the right panel. The regions are ob- 
tained by marginalizing over the channeling fraction. 

For both halo cases, two regions of compatibility are 
present, one for DM masses around 10 GeV and one for 
DM masses close to 50-60 GeV. In the case of heavier 
DM, long-range forces are excluded (the mass bound on 
7770 is set around 100 MeV). In the case of fight DM, 
which is also the overall preferred region, long range 
forces are viable for the whole range of the mediator 
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stands for a threshold of 4 photoelectron and for an isothermal sphere with vq = 220 km s~^ and po ~ 0.3 GeV cm~^. 



masses. However, tlie best agreement is obtained for me- 
diator masses larger than 20 MeV. A 99% C.L. lower 
bound on is about 10 MeV. A variation of the DM 
halo properties does not change dramatically the results, 
as long as the (mean) dispersion velocity is not changed. 
In fact, the comparison between the isotropic (in veloc- 



ity space) isothermal sphere and the anisotropic triaxial 
model shows that the allowed region is slightly enlarged 
in the anisotropic case: both with respect to the DM 
mass and to the mass of the mediator (p. 

In Fig. [lljwe show the bounds obtained in the plane 
e - In this case, we choose to focus on bounds de- 
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FIG. 13: Bounds on the light-mediator mass as a function of the DM mass m^, as obtained by the analysis of the DAMA 
annual-modulation result. The allowed regions are reported for C.L.'s of 90%, 99% and 99.99%. The left panels stands for an 
isothermal halo with velocity dispersion vq = 220 km s"^ and local density po = 0.3 GeV cm"''; the right panels refers to the 
case of a triaxial halo with the Earth located on the major axis [7pJ, with velocity dispersion vo = 220 km s"^ and local density 
po = 0.84 GeV cm"^ 



rived for a kinetic mixing between dark and ordinary pho- 
tons, which is described by a lagrangian term of the form 
f ^^v ^^dLik- These type of coupUngs can produce ex- 
perimentally observable effects in many different physical 
situations. A summary of those constraints are discussed 
e.g. in Refs. [HI], [S2 (and references therein) and are 
reproduced in Fig. 14 as (blue) regions on the left. They 



constrain large couplings and very light mediators. 

Coming to the relevant regions which can be derived 
from the analysis of DM direct detection experiments, 
we first notice that the DM signal depends on four 'dark 
sector' parameters: m^, e, Z', and ajark- Essentially, it 
involves an extra coupling (between DM and dark pho- 
ton, proportional to Z'^adark) with respect to dark pho- 
ton signals in laboratory experiments. Therefore, there 
is more freedom associated, and observational bounds on 
the dark photon properties can only partially constrain 
the allowed parameter space. In this respect, we define 
k = Z' \/ adark/asM and consider two benchmark cases, 
k = 1 and k — 10. In Fig. [T4j we show a combined fit 
to DAMA and CoGeNT annual modulation data. The 
analysis is performed along the line explained in Sec. |III[ 
We combine the two data sets by adding the of the 
two experiments together, and show those regions where 
the absence of signal can be excluded at 8~a C.L.. For 



DAMA, we vary the fraction of channeling in its allowed 
range. From Fig. [l4]we see that in the 'totally symmet- 
ric' case with k — 1 (upper orange region), light media- 
tors are excluded, and only dark photons with > 100 
MeV can simultaneously satisfy the constraints and pro- 
vide a suitable interpretation for DAMA and CoGeNT 
data. However, for k > 10, the whole range of light- 
mediator masses is allowed. It is worth to point out that 
such values of k can be easily predicted in models with 
a strongly coupled dark sector or including 'composite 
DM' with large dark charge. 



Another class of constraints arises from bounds on DM 
self-interaction in cosmic structures, a bound which de- 
pends on the DM coupling with the dark photon and on 
the mass of the dark photon, while it is independent on 
the kinetic mixing parameter e. The long range force be- 
tween DM particles implies that DM is more collisional 
than in the standard WIMP case. A first example of 
such constraints comes from observations of systems of 
colliding galaxy clusters, such as the Bullet cluster |83j . 
which point towards collisionless DM. A robust bound of 
c/wx — -'^•25 cm^/g has been placed on the size of the 
self-scattering of DM [84] . To infer the constraint on the 
mass of the dark photon, we compute the weighted cross 
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FIG. 14: Regions of daxk photon mass m0 and kinetic mixing 
e compatible with DM direct detection experiments, as com- 
pared with existing independent experimental bounds shown 
as (blue) areas on the left (see Refs. |81|.|82) and references 
therein). The upper (orange) and lower (red) regions are 
derived combining DAMA (with the fraction of channeling 
marginalized) and CoGeNT annual modulation data, and re- 
fer to regions where the absence of a signal can be excluded 
with a C.L. of 8-<j. The two regions differ in the choice of the 
value of the effective coupling k of dark matter to the dark 
photon (defined as k = Z' -y/adark/asM): the upper (orange) 
regions refers to fc = 1, the lower (red) region refers to fc = 10. 



section; 



d^vid\2fivi).f{v2) 



,„ da , 



cos 



Q) (27) 



that measures the rate at which energy is transferred 
in the system of colliding clusters. Here j(v) — 



-(v/voY 



is the velocity distribution of dark mat- 



ter particles in the cluster, assumed to be Maxwellian, 
is the velocity dispersion, 9 is the scatter angle in the 
centre of mass frame and the cross section on the right 
hand side refers to the elastic processes. For dark mat- 
ter particles interacting with long range forces, the self- 
interaction is analogous to a screened Coulomb scattering 
with the plasma, which is well fit by the following cross 
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An 

ml 



/32ln(l + /3-i), 
(1 + In /3- 0.5 In"^)' 



13 < 0.1, 
0.1 < /3 < 10^, 
/?> 10^, 



(28) 

where /3 = 2 Z'^ adarlc"T.^/('71xWrol) ^^'^ ^rcl = \vi - W2I 

is the relative velocity of dark matter particles. Consid- 
ering typical velocity of collision in the Bullet cluster of 
4700 km/s, a dark matter mass of 10 GeV (to fit direct 
detection observations) and two values of the parameter 
k ^ (1,10) => Z'^adark = (asM, lOOasM), we get that 
the bound on the self-interaction is exceeded if the mass 
of the dark photon is smaller than (1, 20) MeV. 

A second class of constraints comes from the fact that 
large DM self-interactions causes a rapid energy transfer 
between DM particles and thus tends to drive DM ha- 
los into spherical isothermal configurations [851 - I89] . The 
observation of DM halo ellipticity allows to put further 
constraints on the size of self interaction. To estimate the 
impact that self interaction via screened Coulomb scat- 
tering has in isotropizing the shape of dark matter halo, 
we compute the relaxation time Tj. = F^^^^ where 



d^Vid^V2 f{vi)f{v2) n^Vrci<yivi-ci) {vrd/voY 



(29) 

is the rate at which energy is transferred in the system. 
This relaxation time provides an estimate on the effects 
of self interactions on the dynamics of a virialized astro- 
physical object with number density — p^/m^ and 
velocity dispersion vq- Indeed if it is much longer than 
the age of the object (say r ^ 10^" years for "old" ob- 
jects), we expect that self interaction does not alter the 
dynamics. On the other hand, much shorter relaxation 
times imply that an isothermal configuration tends to 
form and thus such scenario would be excluded by the 
observation of few elliptical halos [551 - I5^ . In the case 
of galaxy clusters (p^ ~ lO^'^ GeV/cm-^ and Vq ~ 1000 
km/s), considering again — 10 GeV and k — (1, 10), 
we get that the relaxation time is always much longer 
than the age of clusters and thus the mass of the dark 
photon is basically unbounded. 

On the other hand, following the same analytic ap- 
proach for virialized halos at galactic scales, where the 
dark matter energy density is higher (p^ ^ 1 GeV/cm^) 
and velocity dispersion is lower [vq — 240 km/s for a 
galaxy and vq ~ 10 km/s for a dwarf galaxy), one can 
in principle put very stringent constraints on m^, ex- 
cluding < 100 MeV. However, Eq. (29 1 is only an 



approximation, and since galaxies form in a fully non- 
linear regime, an N-body numerical simulation explicitly 
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including long-range interactions would be in order to 
properly address the size of the constraint. Moreover, we 
note that interactions with baryons (possibly in the long- 
range regime as well) could substantially affect the halo 
shape. This very interesting subject deserves a dedicate 
treatment which is beyond the goal of this paper. 

V. CONCLUSIONS 

In this paper we have discussed the current positive 
indications of a possible dark matter signal in direct de- 
tection experiments, in terms of a mechanism of inter- 
action between the dark matter particle and the nuclei 
occurring via the exchange of a light mediator, resulting 
in a long-range interaction. We have therefore analyzed 
the annual modulation results observed by the DAMA 
and CoGeNT experiments under the hypothesis of this 
type of long-range interactions, and derived bounds on 
the relevant parameters at hand, namely the DM mass 
and the effective cross section, as well as the mass of the 
light mediator. 

We find that long-range forces are a viable mechanism 
which is able to explain the modulation effects, and we 
have obtained that the preferred range of masses for the 
light mediator is for values larger than a few MeV. This 
is correlated with a DM mass around 10 GeV. We have 
also obtained that the long-range forces mechanism is 
more constrained in the case of large dispersion veloci- 
ties of the DM particles in the galactic halo: in the case 



of dispersion velocities larger than about vq '■^ 250 km 
s~^ the long-range forces are less favored (although not 
excluded). This result is stable over the whole variation 
of the light-mediator mass. We can therefore conclude 
that the annual modulation effect reported by DAMA is 
compatible with long-range forces when the DM mass is 
preferably in the range from 8 GeV to 20 GeV and the 
mediator mass is in excess of about 10 MeV, depending 
on the actual galactic halo model. 
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